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Abstract Gamma rays at rest frame energies as high 
as 90 GeV have been reported from gamma-ray bursts 
(GRBs) by the Fermi Large Area Telescope (LAT). 
There is considerable hope that a confirmed GRB de- 
tection will be possible with the upcoming Cherenkov 
Telescope Array (CTA), which will have a larger ef- 
fective area and better low-energy sensitivity than 
current-generation imaging atmospheric Cherenkov 
telescopes (IACTs) . To estimate the likelihood of such 
a detection, we have developed a phenomenological 
model for GRB emission between 1 GeV and 1 TeV 
that is motivated by the high-energy GRB detections 
of Fermi-L AT , and allows us to extrapolate the statis- 
tics of GRBs seen by lower energy instruments such as 
the Swift-JiAT and BATSE on the Compton Gamma- 
ray Observatory. We show a number of statistics for 
detected GRBs, and describe how the detectability of 
GRBs with CTA could vary based on a number of 
parameters, such as the typical observation delay be- 
tween the burst onset and the start of ground observa- 
tions. We also consider the possibility of using GBM 
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on Fermi as a finder of GRBs for rapid ground follow- 
up. While the uncertainty of GBM localization is prob- 
lematic, the small field-of-view for IACTs can poten- 
tially be overcome by scanning over the GBM error 
region. Overall, our results indicate that CTA should 
be able to detect one GRB every 20 to 30 months 
with our baseline instrument model, assuming consis- 
tently rapid pursuit of GRB alerts, and provided that 
spectral breaks below ^100 GeV are not a common 
feature of the bright GRB population. With a more 
optimistic instrument model, the detection rate can 
be as high as 1 to 2 GRBs per year. 

Keywords gamma rays: bursts • telescopes 



1 Introduction 

The observation of gamma-ray bursts (GRBs) with 
ground-based atmospheric Cherenkov telescopes has 
been a tantalizing possibility in recent years. Pow- 
erful >10-meter telescope arrays such as H.E.S.S., 
MAGIC, and VERITAS have come online in the 
last decade, and satellite detectors such as the Swift 
Burst Alert Telescope (BAT) are capable of pro- 
viding the necessary localization of GRB events 
within seconds over the Gamma-ray burst Coordi- 
nates Network (GC1N0). Despite major campaigns 
to respond to sat ellite burst alerts at all three of 
these instruments ( Aharonian et al 2009 ■ Albert et al 



2007t iGarczarczvk et afcOQStlAcciari et afcoill) . and 
dozens of follow-up attempts, no conclusive detection 
of a GRB with an IACT has yet been made. Air 
shower arrays have also played a complementary role 
in the search for GRBs. A hint of emis sion was de- 
tected by the Milagrito air-shower array ( Atkins et all 
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2003); however no d etections were fou nd by the later 
Milagro experiment ( Abdo et"al 2007 ) . 

Prior to the launch of Fermi on June 11, 2008, 
knowledge about the emission of GRBs above 100 
MeV was limited to a small number of events ob- 
served simultaneously in the EGRET and BATSE in- 
struments on the Comp t on Gamma-ray Obser vatory 
(CGRO) dPingusI Il99a iLe and Dermerl l2009h . One 
fascinating finding by EGRET was the discovery of 
an 18 GeV photon associated with GRB 940217, 1.5 
hours after the event. This was a much longer time 
than the duration of the burst as measured at lower 
energies by BATSE, which determined a T90 of 150 
secondfd (defined as the time between the arrival of 
5 percent and 95 percent of the observed fiuence). 
Though the statistics of these EGRET observations 
were quite limited, they suggested that high energy 
emission in GRBs did occur in some fraction of events, 
and that it could last longer than the lower energy 
emission. 

Fermi -LAT, covering the e nergy range of 20 MeV - 
300 GeV (jAtwood et afeoooh . has now detected emis- 
sion from over 20 GRBs, of the some 80 detected by 
GBM at 8 keV - 40 MeV energy range (jMeeean et all 
20091) . Photons from four of these LAT GRBs were 
detected above 10 GeV, and 2 above 30 GeV. At 
present, the highest energy photons that have been as- 
sociated with any GRB are a 33.4 GeV photon from 
long-duration GRB 090902B, and a 31 GeV photon 
from short-duration GRB 090510. With th e redshift of 
GRB 90902B determined to be z = 1.822 ()Abdo et all 
2009a), this implies a rest-frame energy of 94 GeV. 
The LAT therefore confirms that emission in the 10 
to 100 GeV decade occurs in at least a small fraction 
of both short- and long-duration GRBs. However it 
is not clear how these findings for bright sources ex- 
trapolate to the rest of the population, and whether 
suppression of GeV-scale emission might also happen 
in a substantial number of cases. 

The other major feature of high-energy gamma- 
ray emission seen by the LAT is the verification of a 
timescale for the VHE emission that is often longer 
than that seen by GBM or other experiments sensi- 
tive to soft gamma-rays. An unexpected finding is the 
delayed onset of emission above 100 MeV, t ypically 
by ~ 10 percent of t he GBM T90 duration ([Dermer 
2010). As discussed in iGhisellini et~al (|2010l ). emission 
at high energy is seen to then continue well beyond 
this time, with a lightcurve described by a powerlaw 
with slope -1.5. Understanding the source of this emis- 

2 http : //gammaray .nsstc .nasa. gov/batse/grb/catalog/4b/ 



sion, which begins well within the prompt phase of the 
burst but continues into the afterglow time period, is 
challenging for models of the high-energy production 
mechanism. 

GeV-scale emission could arise through several 
mechanisms, and understanding the impact of each on 
the cumulative spectra will require multiwavelength 
observations over many orders of magnitude in en- 
ergy, combined with high event statistics. High en- 
ergy emission during the prompt phase of the GRB 
can be most simply explained by a spectral exten- 
sion of the internal shock processes (inverse Comp- 
ton ( IC) and synchrotron) that produce the keV-MeV 
flux ( Sari and Piranl 1997 ) . The observed spectrum of 
GRB 080916C, seen over ~ 7 orders of magnitude by 
the Fermi-LAT and GBM instruments, c ould be ex- 
plaine d by a constant synchrotron origin (jAbdo et al 



2009dl ). Other possibilities include emission of GeV 



photons fr om external sho cks in the early afterglow of 
the GRB ( Fan et al 20081 ) or from the reverse shock 
formed whe n the GRB ejecta encounter the interstel- 
lar medium (Wan g et a 12005). The former can explain 
the delayed onset of high-energy emission seen in most 
LAT-detected GRBs. Several authors have preferred a 
purely synchrotron ori gin in the external shock to de- 
scrib e the LAT GRBs (IGao et alll2009l:|Ghirlanda et al 
20101: iKumar and Barniol Dur an 2010), in contrast to 



one in voking inverse C ompton from the shock elec- 
trons ( Zou et al 2009 ). Finally, hadronic processes 
have been proposed as a source of the high en- 
ergy component, an idea that connects this radiation 
to the production of ultra- high energy cosmic rays 
(jRazzaaue et a faiol: lAsano et alll2010i) . 

The redshifts at which GRBs have been de- 
tected span from the local universe to z — 8.2 
(|Salvaterra et all l2009a|) . corresponding to ^95 per- 
cent of the age of the universe. Confirmed redshifts 
for LAT-detected GRBs range span a wide range of 
this distribution from z = 0.736 (GRB 090328) to 
4.35 (GRB 080916C). This suggests that an IACT- 
dctcctcd GRB could occur at essentially any redshift 
where star-formation has been observed. 

The cosmological UV-IR background radi- 
ation produces a barrier to hi gh-energy pho 



tons at extragalactic d istances (jNikishovl [1962 



Madau and Phinnevlll996l ). Moreover, the large ma- 



jority of GRBs are believed to exist at high redshift, 
where the background flux is highly uncertain, and 
gamma-rays at ob served energies as low as 10 GeV 
can be impacted (jGilmore et a l l2009h . The effective 
area, angular resolution, background rejection ca- 
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pabilities of IACTs are strongly energy-dependent, 
and all of these properties decline in quality below a 
few hundred GeV for current-generation instruments. 
Ground-based GRB observations therefore take place 
in an energy regime where both the low-energy 
instrument sensitivity and the impact of cosmological 
background radiation must be carefully taken into 
account if realistic predictions are to be made. 

Our goal in this work is not to comment on the pre- 
ferred emission mechanism for high-energy GRB pho- 
tons, but rather to build a phenomenological model 
that best describes this part of the spectrum, based 
on the limited set of GeV burst detections to date 
and the much larger body of data available from 
lower energy experiments. As the weight of the evi- 
dence in the brightest LAT GRBs does indicate the 
presence of emission mechanisms beyond those pro- 
ducing the prompt flux, we have included allowances 
for both extended temporal components and separate 
spectral components from the Band spectrum in this 
work. In the following sections we show our predic- 
tions for the rate at which CTA will detect gamma- 
ray bursts, using the information available to us from 
Fermi-T AT , lower energy satellite experiments, and 
attempted IACT observations with the current gen- 
eration of instruments. In fJ2J we describe the model 
used in this work. In fJ31 we show results for prop- 
erties of detected GRBs, calculating both the rates 
at which detections occur in our model and the pho- 
ton statistics that can be expected from a detection. 
We use the detection rates by the currently operating 
Swift-BAT and Fermi-GBM instruments as the basis 
for the calculation. We also investigate the effect of 
varying critical input parameters in our model, such as 
CTA energy threshold, background rate, and telescope 
response time delay. In SjS] we show the spectra that 
could be available to CTA from a few sample GRBs 
with various properties. The topic of ^3] is a compari- 
son of our model with the upper limits that have been 
set by observations to date with IACT instruments, 
and with the GRB rate observed with Fermi-h AT . In 
35] we summarize and discuss our findings. 

2 Modeling the Detection of GRBs 

To a large extent, the challenge of modeling gamma- 
ray bursts arises out of the large variance in prop- 
erties seen between events, and the lack of a simple 
model describing the radiative mechanism. In partic- 
ular, each of the 4 bright GRBs seen by LAT above 10 
GeV shows differing behavior. The GeV-scale emission 



in GRB 080916C, observed by Fermi three months af- 
ter launch, was found to be well-described by an ex- 
tension o f the Band functio n seen at keV and MeV 
energies ( Abdo et~al 2009dl ). Separate spectral com- 
ponents from the Band function were found to be re- 
quired to match the GeV-scale emission of the three 
other brightest GRBs in the LAT catalog. Short- 
duration GRB 090510 was found to be dominated 
by a hard spectral component of index -1.62 above 
~ 100 MeV, compar ed to the average Band high- 
energ y index of -2.4 ( Abdo et all 2009bl Ackermann 
2010h . Long GRB 090902B was dominated above 100 



MeV by emission with a spectral index determined by 
LAT to be -1-93. comp ared to an upper Band index of 
-3.8 (jAbdo et alll2009dY The high energy emission also 
extended in time well past the prompt phase as deter- 
mined at lower energies, with an only slightly softer 
spectrum of -2.1 on a timescale of 10 3 s. Finally, GRB 
090926A, detected up to 19.6 GeV with the LAT, was 
best fit with a high energy hard component of spectral 
index of -1.72 and an exponential cuto ff at 1.4 GeV, 
with a high-energy Band index of -2.63 ( Breeeon et all 
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An earlier work, iGilmore et al ri2010lh addressed 
the question of detection of GRB photons with the 
MAGIC telescope and the Fermi-LAT and the impact 
of the UV-optical background light. The basis for this 
model was the population of GRBs seen by Swift-BAT, 
for which redshifts have been confidently determined, 
and the high-energy statistics of GRBs detected with 
CGi? 0-EGRET. The flux at the MAGIC energy range 
was calculated by assuming a power-law spectrum 
continuing to ~ 200 GeV. Due to the absorption by 
the extragalactic background light (EBL), photons at 
higher energy were greatly attenuated and insignifi- 
cant in number, and the UV-optical EBL therefore has 
a large impact on GRB dctectability. This calculation 
showed that MAGIC was capable of detecting tens to 
thousands of photons from a bright GRB at high red- 
shift, provided that the burst could be observed with 
a sufficiently low energy threshold. 

In this work, we attempt to improve on this pre- 
vious calculation to make predictions for CTA. To 
construct a model for the distribution of GRB prop- 
erties, we draw upon keV and MeV data from two 
satellite experiments: CGRO-BATSE and Swift-BAT, 
plus data from the high-energy GRB detections with 
Fermi-LAT. While Swift has provided us with a large 
number of GRBs with determined redshifts, the BAT 
instrument is generally not capable of resolving the 
Band function peak of the GRB spectra. We there- 
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fore combine the Swift redshift distribution with the 
distribution of Band function parameters seen by the 
BATSE experiment on the CGRO satellite, at ob- 
served energies between 20 keV and 2 MeVQ Band 
function fits for BATSE GRBs are taken from the 
BATSE 5B catalog (Goldstein et al. in prep). To ex- 
trapolate the spectra of these bursts to VHE energies, 
we use the statistics of GRBs seen by the Fermi-LAT 
instrument above 100 MeV within its first two years 
of observations. 



2.1 High energy extrapolation 

Predicting the GeV-scale emission of GRBs from the 
well-sampled statistics of lower-energy instruments re- 
quires a considerable amount of extrapolation. Some 
4 logarithmic decades in energy lie between the upper 
extent of the BATSE energy range (~ 2 MeV) and 
the energy threshold of CTA, which we consider to 
be between 10 and 25 GeV. We describe our two dif- 
ferent approaches to performing this extrapolation in 
this section. 

2.1.1 Band-function extension model 

As a minimal model, we consider the fluence predicted 
for GRBs at high energy without any sign ificant devi- 
ation from the Band fit (jBand et all ll993): 



— oc E a e~ E( - 2+a)/E » 
dE 

dE y 2 + a pJ 



a-B 
E < - — -E„ 



a-P e l3-a. 



a-B 
E > - — -E„. 



(1) 



Here a and /3 are low- and high-energy spectral in- 
dices, and E p is the "peak energy" describing the lo- 
cation of the turnover. In this extended Band-function 
model, termed "bandex" in subsequent plots and dis- 
cussions, the high energy spectrum is assumed to 
merge seamlessly with the spectral fit determined at 
lower energy. A similar model was used in estimat- 
ing the de t ectab ility of GRBs with Fermi LAT in 
Band et all (|2009l ). The high energy normalization is 
therefore determined by the Band function peak en- 
ergy and normalization and the upper energy index 8, 
which continues to GeV energies. In Figure[TJ we show 
the distribution of values for 8 against the BATSE 
fluence. A minority of the GRBs in the sample, about 

3 We preferred using the BATSE catalog instead of the 
GBM catalog because of its better instrument sensitivity 
and the much larger number of bursts detected. 

4 http : //gammaray .nsstc .nasa. gov/~goldstein/ 
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Fig. 1 Band function high-energy indices (0) vs BATSE 
fluence for the GRB population, along with distributions 
for each parameter. 



13 percent, have a hard spectrum with 8 > —2. We 
have enforced the requirement that the total fluence 
per logarithmic decade not be higher in the GeV range 
than at BATSE energies, and thus we reset these cases 
to have 8 = — 2. This requirement is consistent with 
the LAT-GBM fluence relation for long GRBs ob- 
served by LAT, but not the short bursts (see Figure 
(2), where fluence ratios greater than 1 have been ob- 
served. 

2.1.2 Fixed parameter model 

In this fixed-parameter ( "fixed" ) model, we make the 
assumption that the relative fluence between BATSE 
energies and GeV energies can be described by a single 
ratio, which we set here to 0.1. The choice of this pa- 
rameter is based upon the corresponding ratios found 
for simultaneously-observed BATSE-EGRET GRBs 
and for GBM-LAT long -duration GRB s. In Figure [2j 
we reproduce figure 1 in Dermer ( 2010h . showing this 
relation for several Fermi and CGRO GRBs. The spec- 
tral index at high energies is set to —2, consistent with 
the mean value for EGRET GRBs of 



1.95 (jDingus 

19951 : iLe and Dermerll2009h . and near the center of the 
distri bution for LAT-detected events ( Ghisellini et all 
2010). In general this model requires a significant de- 
parture from the extrapolated Band function, and im- 
plies the appearance of a separate high-energy spec- 
tral component. As discussed above, such components 
were seen in GRB 090902B and GRB 090510. A sepa- 
rate spectral component was also the preferred model 
in describing the total time-integrated emission from 
GRB 090926A, albeit with a spectral cuto ff of this 
component at 1.4 GeV ( Breeeon et all 2011 ). Though 
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Fig. 2 A comparison of the LAT and GBM fiuences 
observed in L AT-detected GRBs, reproduced here from 
iDermed (|20ld 1. Black squares are the result for long- 
duration GRBs, while red diamonds are short GRBs. 
Dashed lines indicate LAT-GBM fluence ratios of 0.1, 1.0, 
and 10.0 (bottom to top). Stars are the analogous EGRET- 
BATSE fluence relations for bursts detected by EGRET on 
CGRO. 



this cutoff component for GRB 090926A was found 
with high significance in the integrated fluence, in 
time-resolved analysis this fit was only preferred over 
simpler power laws within a single narrow time win- 
dow. As GRB 090902B was found to have a fluence 
ratio of nearly 10 percent and a time-integrated h igh 
energy index 1.93 (|Abdo et afe009at iDerme rl l2010f) . it 
can be considered the prototypical GRB in motivating 
this extrapolation scheme. While the GeV-BATSE ra- 
tio observed for short GRB 090510 was considerably 
higher, ~ 1, we have not included a separate account 
of the short population because the emission of these 
GRBs at multi-GeV energies remains very poorly un- 
derstood, and they are a small part of the Swift GRB 
sample, ~ 9 percent. Additionally, the time delay to 
see GRBs from the ground makes short GRBs very 
difficult to detect even with a flux factor of unity, and 
including such an possibility is found to have little 
effect on our findings. 



2.1.3 High-energy lightcurve 

As we shall see, the lightcurve and emission duration 
at GeV energies are critical variables in determining 
the detectability of GRBs from the ground, where the 
response time of the telescope to transient alerts from 
satellite instruments limits observations to Si 1 min 
after the start of the event, and the background can 
obscure low-luminosity emission occurring over long 



timesc ales. Motivated by the finding of Ghisellini et all 
( 2010l ). we will assume that the GRB lightcurve in the 
early afterglow phase can be described as a power- 
law falloff. The prompt phase of the GRB can be 
demarcated by T90 as determined by lower energy 
gamma-ray instruments. Luminosity during this phase 
is often seen to fluctuate rapidly and unpredictably, 
with spiked emission features tha t undergo rapid ex- 
ponential time decay ( Piran [20041 ) . Since only consid- 
ering time-averaged behavior (where the typical er- 
ratic prompt emission of GRBs is neglected) will not 
affect our results in term of detection rate and pho- 
ton statistics, we describe this phase as having con- 
stant flux. Our total modeled GeV lightcurve then is 
a plateau from the burst onset (t = 0) to t = T90, 
followed by a powerlaw falloff: 



F(t) = F ; 



t < T90 (2) 
t > T90 



Here 7 is the power-law index of the afterglow 
lightcurve. We will use 1.5 as a fiducial value, but 
will also explore the impact of other possibilities in 
34.1.21 Note that for this value, two-thirds of the total 
emission emerges after T90. We assume no spectral 
evolution between the prompt and afterglow phases. 



2.2 Redshift distribution 

The observed fluence distribution for GRBs is not as- 
sumed to be directly dependent on redshift. However, 
redshift is a crucial factor in determining GRB de- 
tectability because the cosmological opacity due to 
EBL is determined by the GRB redshift. We make the 
assumption in this calculation that the redshift distri- 
bution of GRBs to which CTA responds will be simi- 
lar to that seen by the Swift-3 AT experiment, which 
is the only large sample of GRB redshifts available. 
Approximately one-third of the GRBs in the Swift 
population have well-determined redshifts. In Figure 
[21 we show the distribution of Swift redshifts for 167 
GRBs, taken from the online Swift GRB Lookup Ta- 
blc0, along with our fit to the distribution which is 
used in most of this analysis. In £14.21 we will make a 
speculative alteration to this distribution to describe 
the redshift distribution of GBM GRBs. 

In this work, we make the assumption that red- 
shifts and observed fluence are uncorrelated. In Fig- 
ure 21 we show the distribution in fluence for Swift 

5 http: //swift .gsf c .nasa. gov/docs/swif t/archive/grb_t able/ 
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Fig. 3 The redshift distribution determined for Swift GRBs. 
The dashed line shows the fit used in this work. 



GRBs divided into tertilcs in redshift. The distribu- 
tion is not found to evolve strongly in redshift, and 
the lowest redshift bin actually h as the lowest median 
fluen c e. It has been suggest ed ( Llovd-Ronning et~al 



2002; ISalvaterra et all l2009bl ) that luminosity evolu- 
tion in redshift (e.g., by a factor (1 + z) a with a £ 1) 
is required to best fit the redshift-luminosity relation 
seen in GRBs. While the existence and possible origins 
of such a factor remain controversial, such an evolu- 
tionary term could account for our findings in Figure 

H 



2.2.1 CGRO-Swift fluence matching 

The Swift-BAT population of GRBs is found to have 
a lower average fluence distribution than the Band- 
function fits CGRO-BATSE 5B population that is 
sampled to determine GRB fluence and spectral prop- 
erties for our burst samples, when the latter is inte- 
grated over the BAT energy range. We have adjusted 
the total fluence of the BATSE burst population by 
a global factor to better match that of Swift, which 
we consider as the GRB trigger instrument in 34.11 
and in 34.3.11 by using a Kolmogorov-Smirnov test to 
minimize the difference between the distributions for 
the brightest 50 percent of bursts in the BATSE sam- 
ple and the brightest 50 percent in the BAT sample. 
The fact that we have restricted our fit to the bright- 
est 50 percent of GRBs is motivated by our finding 
that GRBs with less than median fluence are not de- 
tectable by CTA even under optimal conditions. The 
multiplier applied to BATSE fluences in 34.1l is 0.501. 
The distribution of the adjusted BATSE population 
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Fig. 4 The cumulative fluence distribution of Swift GRBs in 
three redshift bins, which have been chosen to each contain 
an approximately equal number of events. Blue line: z < 1.2, 
green line: 1.2 < z < 2.7, red line: 2.7 < z. The dotted line is 
the integrated fluence in the total BATSE population over 
the BAT energy range, 15 to 150 keV, after the correction 
discussed in £|2.2.1I Note that this fluence distribution is for 
all GRBs, not only those with redshift, and therefore it is 
not necessarily expected to match the Swift distributions as 
GRBs with redshift are found to be slightly brighter than 
the population as a whole (Fig. [5]). 



after this correction is shown alongside the Swift-BAT 
fluence distributions in Figure [H 



2.2.2 EBL attenuation 

The EBL, specifically at UV-optical wavelengths, is 
responsible for attenuating the signal of high energy 
gamma-rays. In some EBL models, this attenuation 
can affect gamma rays at observed energies as low 
as 20 to 25 GeV for high redshift sources. The ef- 
fect of the EBL is to reduce the number of gamma- 
rays received at high energy, and to reduce the de- 
tectability of high-redshift GRBs. We use as a stan- 
dard assu mption in this work the EBL model and 
opacities of lGilmore et all (1201 lh (GSPD11), based on 
the semi-analytic modeling of ISomerville et all (|201ll ) 



and Somerville et all 



(2008). However, large uncertain- 
ties in the EBL normalization are unavoidable at high 
redshift, and it is useful to see exactly how this uncer- 
tainty can influence our predictions. In 34.1.31 we will 
look at our results assuming a few different models for 
the background light. 



2.3 Telescope properties 

The second step in constructing our model is a 
parametrization of the performance of the CTA. As 
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Fig. 5 Integral distibution of Swift-observed GRBs, com- 
pared with the BATSE population integrated over the same 
energy range. The red dashed line shows the BATSE popu- 
lation from 20 keV - 2 MeV, and the solid curve shows these 
same GRBs integrated over the Swift-BAT energy coverage 
of 15 - 150 keV. The dotted red curve shows the BAT-band 
distribution of the BATSE GRBs after the proposed ad- 
justment to best match the distribution of the brightest 50 
percent of Swift-BAT GRBs (solid blue). Solid and dashed 
blue lines are the Swift-BAT fluences for all GRBs and those 
with redshifts, respectively. 



many of the array properties are indeterminate at the 
time of writing, we have r elied on the design concept 
for the array described in lActis et all ([20111 ) , as well 
as reasonable extrapolations from the current gener- 
ation of IACTs, particularly the MAGIC and VERI- 
TAS telescopes. 



2.3.1 Effective area 

Our assumptions about the effective area of CTA 
are based on Configuration E, which assumes a cen- 
tral cluster of four 24-meter class large-size telescopes 
(LSTs) that provide sensitivity to the lowest energy 
gamma-rays, and an additional 23 medium-size tele- 
scopes of the 12-mctcr class (MSTs) providing sensi- 
tivity at higher energies, £ 100 GcV. Sensitivity at 
energies above 1 TeV, which is provided by more dis- 
persed arrays of 7-meter class small-size instruments 
(SSTs), is not crucial to our results here, as most 
GRBs will occur at redshifts for which emission at 
these energies is strongly attenuated by the EBL. 

The effective area function of the instrument, after 
all analysis cuts have been performed, is used to deter- 
mine the counts per GRB and the significance of the 
detection presented in the next section. As the actual 
function is unknown at this point in time, we assume 
two functions for the effective area, which are shown 



in Figure El Each of these functions includes contri- 
butions from the LST and MST arrays, which dom- 
inate the total effective area below and above ~ 100 
GeV, respectively. For the more conservative of the 
two, labeled "CTA baseline" , the LST contribution is 
created by shifting the standard VERITAS area func- 
tion to the 25 G eV threshold that is expected for CTA 
(|Teshimall201lT) . The MSTs use a function that is un- 
modified in the energy dimension, but has a normal- 
ization factor that assumes a linear scaling in effective 
area with telescope number; we adopt a factor 25/4 
= 6.25 for this case. We believe these numbers to be 
a reasonable estimate of the capabilities of the array. 
As an alternative, we also present results that incor- 
porate several enhancements to the baseline assump- 
tion, and are intended to represent the best possible 
performance that can reasonably be expected from the 
instrument. For this case, labeled "CTA optimistic", 
the LSTs have the same normalization as the CTA 
baseline array, but with an energy threshold of 10 
GeV, which might be achieved through improved trig- 
ger and background rejection techniques. The MSTs 
are given a normalization 3 times that of the baseline 
(75/4 = 18.75). This could either be taken to repre- 
sent the coverage of a 75 telescope array, or a smaller 
number of telescopes if the effective area increase with 
telescope number is found to scale at a faster than 
linear rate. The energy shifts for the LSTs in each 
case are assumed to take place multiplicatively, i.e., 
A C J S A {E) = AY/f R (k * E), with k = 4 for the baseline 
area function and 10 for the optimistic. 

These parametrizations of the telescope effective 
area refer to a source at zenith. The changes in tele- 
scope performance away from zenith, i.e. at an angle 
9 Z , are considered using adjustments to the effective 
area function. Viewing sources at increasing 6 Z gen- 
erally introduces a higher energy threshold to the ob- 
servation, due to the increasing distance between the 
shower core and the telescope array. The increased 
area of the light pool is also the reason why the effec- 
tive area increases at higher energies. The first effect 
is parametrized as a multiplicative energy shift in the 
effective area function; this is considered as a shift in 
the telescope energy threshold by a factor cos(6 z )~ 3 
away from zenith. The choice of this index, and the use 
of a multiplicative shift rather than an additive shift, 
are both motivated by a fit to the VERITAS effective 
area at various elevations, with cuts optimized to soft 
sources (VERITAS Collaboration, private communi- 
cation). The second effect is treated by enhancing the 
effective area at all energies by a multiplier cos(# z )~ 2 , 
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1000 



Fig. 6 The effective area functions used in this work. Solid 
red is the VERITAS effective area with standard cu ts , and 
the dotted blue line is the MAGIC i jAlbert et all2008t ) imple- 
mentation with standard cuts, shown here for comparison. 
The two dotted black curves are the effective area functions 
used in this work, denoted CTA realistic (lower) and CTA 
optimistic (upper). 



which represents the geometrical increase in the area 
of the Cherenkov light pool. The off-zenith effective 
area is then 



A?L A (E,e z ) = ^A c J s A {^E,e z = o), 



x eff 



(3) 



where r\ — cos(# z ) and A^^ A {E, 9 Z — 0) are the at- 
zenith functions shown in Fig. [51 

2.3.2 Instrument background 

Understanding the instrument backgrounds that will 
impact GRB observations is critical to predicting the 
detection rate. Unfortunately, while published rates 
from present-day telescopes can be used to predict the 
background that will impact observations at ^100 
GeV, little is known about how these rates will extend 
to lower energies. 

We base the background rate in our analysis on 
that of typical VERITAS observations, using the as- 
sumption that the four LSTs will achieve a similar 
rejection fraction as the four VERITAS instruments 
above each of their respective energy thresholds, after 
assuming a power law scaling of the background from 
100 GeV to lower energies. Meanwhile, it is assumed 
that the MSTs will, at a maximum, have the same 
background as the VERITAS array multiplied by the 
scaling in effective area normalization. The MSTs con- 
tribute only a minority of the total background rate in 
our analyses, so an overestimate of their contribution 
would have little effect on our results. The background 



spectrum is assumed to be a power law of index -2.7. 
We then take the background at a given energy to 
scale in proportion to the effective area at that energy. 
The total background rate for a telescope set (LST or 
MST) is then 



dN\ 
dt) bkg 



A 



E Aejf(E) dE. 



(4) 



The normalizing factor A is chosen so as to produce 
the rate for the VERITAS instrument, when this in- 
tegral is performed over the VERITAS effective area 
function. A rate of 6 counts per minute (0.1 Hz) is 
assumed for VERITAS in the case of the baseline ef- 
fective area, and this is reduced to 2 counts per minute 
for the optimistic case. 

Using the baseline effective area curve and back- 
ground rate, together with the scaling described, we 
find a differential sensitivity above 100 GeV similar 
to that sh o wn fo r configuration E in Figure 24 of 
Actis et all (|201lh . Below 100 GeV, we find that as- 
suming a ~35 GeV threshold for the LST array in 
our model produces differential sensitivity similar to 
config uration E , while the 25 GeV threshold we adopt 
based iTeshimal (|201l[ ) gives a sensitivity as much as 
three times better at some energies. 

2.3.3 Response time 

The transient and random nature of GRB emission 
represents the main difficulty in detecting emission 
from these sources. The onboard satellite localization 
time of the event, transmission of the data to the 
ground, the observer's response time, and slew time 
for the IACT all contribute to a total delay time for 
the commencement of observation, which we quan- 
tify in this work as T^eiay The localization time is 
dependent on the instrument and brightness of the 
GRB, but times of < 15 sec are typical. The trans- 
mission time of GRB coordinates is exp ected to be 
nearly instantaneous (jBastieri et alll2005l ). The LSTs, 
which provide coverage at the crucial low energies, are 
expected to have a slew time of 20 to 30 seconds, while 
the MSTs may be somewhat slower. As a standard as- 
sumption, we assume a total response time of 60 sec- 
onds in this work for the LSTs and 100 seconds for the 
MSTs, but will also discuss in the next section the ef- 
fect of varying this parameter. To date, most observa- 
tions with the MAGIC telescopes have commenced af- 
ter considerably longer times despite the instrument's 
rapid slew capabilities, with only a minority occurring 



with total delay times of < 100 sec (jGarczarczvk et a 



GRB prospects for CTA 



9 



2008). It may be that the longer delay times are due 
to reasons other than the mechanical capabilities of 
the instrument. While the inner telescopes of CTA 
are expected to have generally the same slewing ca- 
pabilities as those of MAGIC, we will allow for the 
possibility that future improvements to the GCN and 
telescope alert procedures and observer response time 
could lower the typical delay time from current values. 



3 Limits from current experiments 

Before discussing our results for simulated CTA ob- 
servations, it is useful to compare the predictions of 
our models to the findings of current-generation GeV- 
scale experiments. This enables a test of our model 
beyond the bright events seen by the Fermi-LAT . 



3.1 Analysis of VERITAS GRBs 

While upper limits on GRBs have been published 
for all of the major current-generation IACT exper- 
iments, a problem arising when comparing to these 
limits is the dependence on a particular set of assump- 
tions about the high-energy spectrum, lightcurve, and 
EBL model. These must be equivalent for a mean- 
ingful comparison to be made between any two sets 
of observations and/or predictions. A confirmed GRB 
redshift is also necessary to determine the impact of 
EBL on attenuated flux, which can change flux pre- 
dictions by orders of magnitude. 

Here we consider the GRB li mits from VERITAS 
presented in Acciari et all ( 2011 ). This work consid- 
ered 16 GRBs observed by VERITAS over a 30-month 
period, 9 of which have spectroscopically-confirmed 
redshifts. The analysis assumed a characteristic af- 
terglow decay of to find the optimal integration 
timcscale for detection, as we have applied throughout 
this work. The eff ect of the EBL was c ompensated by 
using the model of Gil more et a 3 (l201lh . as used in this 
paper. Results are presented for gamma-ray spectra of 
dN/dE ~ E- 2 - 5 and E~ 3 - 5 . 

To compare our predictions with these results, we 
have performed a calculatio n of the high e nergy emis- 
sion for the 7 GRBs from lAcciari et all (|201lf) that 
have both redshift determinations and have been ana- 
lyzed assuming a t~ 15 afterglow decay. To model these 
GRBs in such a way that a direct comparison is pos- 
sible, we use a modified version of our fixed model. 
Unfortunately, the flux information provided by Swifb- 
BAT is generally not sufficient to resolve the spectral 



peak of the GRB emission and determine the Band 
function spectral fit. In our fixed model, we assume 
that the flux seen in the Swift-BAT bandpass (15- 
150 keV) is related to the flux in the 20 keV -2 MeV 
band by applying th e common Band function used in 
Gilmore etaf (|2010f ). which leads to a ratio ~ 5 be- 
tween flux in the 20 keV -2 MeV and Swift bands. 
While this facto r is intermediate to the range seen in 
BATSE GRBs (|Preece et allbOQOl) . the considerable 
variations present in the spectral indices and peak 
energy could change this ratio by a large factor. To 
match the spectral index of -2.5 used in the standard 
analysis, we assume a spectral index of -2, as used 
in our work, with a spectral turnover to -2.5 at the 
energy threshold of each GRB observation. 

In Table [TJ we compare predictions from our model 
for GRBs with upper limits placed by VERITAS. In 
one instance, GRB 070521, highlighted in the right- 
most column, the predicted flux exceeds the upper 
limit set by VERITAS by a factor of about 1.5. We 
do not believe that this case alone poses a problem for 
our fixed model, as this factor is much smaller than 
the scatter seen in MeV-GeV fluence ratios for bright 
GRBs (Fig. [5]), and there is an additional degree of un- 
certainty in extrapolating the Swift-BAT bandwidth 
to the BATSE energy range. We also note that the 
Tdeiay and optimal T b s reported in this case are both 
quite high, 1118 s, and 1809 s, compared to the T90 
duration of 38 s. The predicted high energy emission 
for this GRB in our model is therefore reliant on the 
extended-duration lightcurve, up to a timescale of ~ 1 
hour. Finally, in r eporting the redshi ft for GRB 070521 
listed in Table [Q Irlattori et al (|2007l) noted that the 
detection could be spurious due to a faint afterglow 
from the supposed host. If the GRB were at higher 
redshift, then our predicted flux would be lower and 
the disagreement lessened or removed entirely. 



3.2 GRB detection with Fermi-LAT 

Fermi-LAT provides the largest set of high-energy 
GRB detections, and we have compared the rate of 
detections with this instrument with those predicted 
in our models. To describe th e LAT, we ha ve used 
the P6 parameters described in Rando (2009) and on 
the Fermi-LAT performance websit^E Our assumed 
effective area is based on the "transient" event type, 
and we assume a background rate of 0.05 Hz within 



6 http: //www-glast . slac . Stanford. edu/sof tware/IS/ 
glast_lat_perf ormance . htm 
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GRB ID 


F BAT 


Redshift 


VERITAS <F> UL 


< F >fixcd 


< F >fixcd/< F >UL 


070419A 


5.58 


0.97 


2.3 x 10- 11 


9.36 x 10- 16 


4.07 x 10" 5 


070521 


80.10 


0.553 


7.6 x 10~ 12 


1.13 x lO" 11 


1.48 


080310 


23 


2.43 


2.8 x 10" 11 


3.27 x 10- 17 


1.17 x 10- e 


080330 


3.4 


1.51 


3.8 x 10" 11 


2.16 x lO" 14 


5.69 x 10~ 4 


080604 


8.0 


1.416 


3.1 x 10" 11 


8.76 x 10- la 


0.028 


080607 


240 


3.036 


9.8 x 10" 11 


5.50 x 10- le 


5.61 x 10~ H 


090418A 


46 


1.608 


6.9 x 10" 11 


5.31 x 10- 12 


0.077 



Table 1 Comparison of VHE flux computed for 7 VERITAS GRBs using a modified version of our fixed model (see text) 
with the upper limits from observation. Columns show the GRB number, the Sure/f-BAT fluence in units of 10 -7 erg cm -2 
s _1 , and redshift. The last 3 columns are VERITAS upper limits, predicted fixed model flux averaged over the observation 
time, and the ratio of the two. Cases where predicted flux exceeds the upper limit are shown in bold. 



the point-spread function of the instrument. Both the 
background rate and the effective area at all ener- 
gies are assumed to evolve function of boresight 
angle uniformly in proportion to the background at 
10 GeV. This analysis uses the same parameters as- 
sumed previously, the only change being that there is 
no delay time in LAT observations, and T a b s (the time 
over which the signal is integrated) is considered on 
timescales as short as 0.1 s. 

We find detection efficiencies (fraction of GRBs 
that are detected) above 100 MeV of 12.6 and 5.1 
percent for the bandex and fixed models, respectively, 
for GBM bursts occurring within 70 degrees of the 
LAT boresight. The specifics of the calculation used 
to find detection efficiency will be described in detail in 
< j4.1l The fact that these results are inverted from the 
pattern seen in CTA results in the following sections, 
where the fixed model generally has a higher detection 
efficiency, can be explained by the lower energy range 
covered by the LAT, which favors the softer GRBs in 
the bandex sample. Th ese fractions ca n be compared 
to the 2-year results of iBissaldil (|201ll) . which report 
270 GBM GRBs within this angle, and 18 LAT de- 
tections; an overall ratio of 6.7 percent. We have not 
accounted for autonomous repoints done by the tele- 
scope, which have occurred in 45 cases, and could po- 
tentially have the effect of raising the detection rate. 
GRBs occurring at the center of the LAT field of view 
are found to have a detection efficiency of about 1.3 
times that of all GRBs within 70 degrees of boresight. 
As spacecraft repointings have only occurred in a rel- 
atively small fraction of cases (45 out of ~540 GRBs) , 
we conclude that the overall impact of these repoints 
on the detection efficiency is expected to be minor, 
even before the observational time delay introduced 
by the telescope slew is taken into consideration. 

As discussed in the Introduction, 4 of these de- 
tected GRBs have had detected emission above 10 



GeV. In both of our models we find that about 30 
percent of LAT-detected GRBs have at least 1 de- 
tected photon above 10 GeV, within 1000 sec of the 
event onset. This is only slightly higher than the cor- 
responding 2-year ratio of 4/18 « 22 percent. Due to 
the long timescale assumed, the model may overesti- 
mate the detection rate in some cases due to practical 
observing constraints, such as instances in which the 
GRB happens to vanish below the horizon. 

4 Results 

The GRB detection capabilities of CTA can be de- 
scribed as the product of two independent factors: 



Detection Rate = DE x TR. 



(5) 



Here DE denotes the detection efficiency, or proba- 
bility that a randomly-selected GRB for which CTA 
is able to take data will be detected with a signif- 
icance of more than 5 standard deviations, and TR 
is the trigger rate at which the telescope is able to 
successfully respond to triggers from satellite instru- 
ments. The product of the two is the rate at which 
confirmed, statistically-significant detections of GRBs 
will take place. The factor TR can be decomposed into 
several independent parameters, which are addressed 
in g~3 

In the following two sections, we show the results of 
our modeling of the detection efficiency, and its depen- 
dence on various instrumental properties. In i j4.ll we 
will assume GRBs follow the statistics seen in Swift de- 
tections, and use the redshift distribution and CGRO- 
Swift flux multiplier of 0.501 that was motivated in 
[2J In i )4.2[ we will consider the alerts provided by the 
GBM instrument on Fermi, and we will update the 
redshift distribution and fluence multiplier to better 
fit this data set. A particular difficulty arising from 
GBM alerts is the large positional uncertainty, which 
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Fig. 7 Basic statistics of GRBs observed with CTA. In each panel the solid black line is the result for the direct extrapolation 
of Band functions (bandex model), from the distribution seen in BATSE GRBs, and the broken blue line is for the fixed 
model, using parameters described in ^2.1.11 and 12. 1.2| respectively. Thick lines are for the baseline effective area function, 
thinner lines are the optimistic function. Left: The integral distribution of sigma values (significance of the source counts) 
for all simulated GRBs in the population. Right: The integral distribution of source photon counts for GRBs which are 
detected, in the timescale bin with maximum sigma. 



is in many cases larger than the CTA field of view, 
and we devote most of this section to addressing some 
possible strategies to maximize the usefulness of GBM 
alerts. 



timescale. The GRB is then assumed to be detected if 
the significance is more than 5 sigma. 



.1.1 Distributions 



id AL 



4.1 Simulated observations and detection efficiency 
for Swift-like GRBs 

For each of the two spectral models considered in W2.ll 
we consider observations of GRBs randomly placed 
within a disk of 75° around zenith. Once the spectrum, 
lightcurve, and telescope effective area and thresh- 
old energy have been determined, we calculate the 
total integrated counts and background counts over 
17 observation timescales (T j, s ), with equal logarith- 
mic spacing from 1 to 10 4 seconds. Observations with 
the telescope are assumed to commence at a time 
T delay after the beginning of the burst, and end at 
Tdeiay + T oos . Observed energies are considered from 
1 GeV to 1 TeV, although most GRBs experience a 
spectral cutoff at energies lower than 1 TeV due to the 
EBL. For each timescale, we calculate the significance 
cr of the GRB det ection, using the m ethod described in 
equation (17) of Li and Ma ( 19831 ). For the purposes 
of this analysis, an on target - off target time ratio of 
1/3 is assumed. The calculated values for a for each 
timescale are compared, and the highest significance 
for the bins that have more than 10 photon counts is 
chosen as the significance for detection of the GRB, 
and the corresponding T f, s is designated the optimal 



In Figure [7] we show the basic statistical results of sim- 
ulated CTA observations for a calculation using the 
effective area curves of Figure [51 Tdeiay = 60 seconds, 
and a maximum angle from zenith of 75°. A majority 
of observed GRBs in both models (~ 90 percent in the 
bandex model; ~ 80 percent in the fixed) do not lead 
to a signal of any appreciable significance (a < 1) for 
the baseline effective area. The detection efficiency in 
the bandex and fixed models is found to be 7.3 percent 
and 11.4 percent, respectively, for the baseline effec- 
tive area, and 16 and 33 percent for the optimistic 
effective area. Overall, the bandex model shows a flat- 
ter distribution of a- values than the fixed model. This 
is due to the additional degree of freedom introduced 
in this model by considering the upper Band index 
in determining the high-energy GRB output in addi- 
tion to the BATSE fluence, leading to a wider range of 
values for the overall high-energy normalization. For 
the same reason, while the bandex model is more pes- 
simistic in its predictions for detection efficiency, the 
detected GRBs in this model do often produce more 
photon counts than the fixed model. This can be seen 
in the right-hand panel of Figure which shows the 
distribution of photon counts for detected GRBs. 

Fig.[5]shows the integral distribution for both total 
source photon count and integral counts above several 
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Source Counts (N) 

Fig. 8 A breakdown of the number of photons seen in de- 
tected GRBs, as a function of energy threshold. This result 
is for the bandex model. The solid curve is the integral dis- 
tribution for the number of all photons. From top to bot- 
tom, dashed lines shown the distribution of the number of 
photons above 50, 100, 200, 300, and 500 GeV. 



Instrument 


DE (bandex) 


DE (fixed) 


CTA (baseline) 


0.0744 


0.115 


CTA (optimistic) 


0.163 


0.328 


CTA (baseline; LST only) 


0.0732 


0.110 


CTA (baseline; MST only) 


0.0231 


0.0310 


VERITAS (E th = 65 GeV) 


0.0241 


0.0281 


VERITAS (E th = 100 GeV) 


0.0216 


0.0235 



Table 2 Summary of detection efficiencies for several in- 
strumental arrangements. In the 'LST only' and 'MST only', 
the effective area and background contributions of the MST 
and LST components are respectively set to zero. We also 
show results for the VERITAS effective area, assuming two 
different energy thresholds. 

energy thresholds for the bandex model with the base- 
line effective area function. The majority of photons 
for most GRBs are seen to arrive below 100 GeV, with 
a significant fraction below 50 GeV, despite the much 
larger effective area provided by the MSTs at higher 
energies. Also, for a majority of detected GRBs the 
expected number of source counts above 300 GeV is 
less than 1. 

Table [5] summarizes the detection efficiencies found 
for a variety of different possibilities in instrument 
configuration and assumed GRB population model. 
In £14.31 we discuss how these results for detection ef- 
ficiency (fraction of GRBs viewed by the instrument 
that will be detected) can be converted into a detec- 
tion rate. We also show results for the VERITAS ar- 
ray in the table, for comparison. In this case we follow 
the same analysis procedure as for CTA, using the 
VERITAS effective area function shown in Fig. El and 
assuming a delay time of 100 s, rather than 60 s. 



Figure [5] shows the expected distribution of red- 
shifts for detected GRBs, compared to the whole pop- 
ulation. The CTA effective area, with sensitivity be- 
low 50 GeV, potentially allows GRB detections at high 
redshift, though those at lower redshift will generally 
have better photon statistics and will therefore be fa- 
vored. Assuming the baseline effective area, few GRBs 
are detected above redshift 2, due to the strong impact 
at higher redshift of the UV- optical EBL at energies 
above 50 GeV. When the optimistic effective area is 
assumed, a subset of GRBs (~ 0.1) are bright enough 
from 10 to 50 GeV to be detectable even at very high 
redshift. These detections are still a minority of the 
full set of detected GRBs, however, and are entirely 
dependent on the low energy performance of the LST 
array. In all cases, the distribution is significantly bi- 
ased towards lower redshifts relative to the Swift dis- 
tribution as a whole, with median redshifts of z = 0.9 
and 1.2 for the baseline and optimistic effective area 
functions, respectively. 

In Appendix [Al we show a number of other prop- 
erties for GRBs that pass the detection criteria. This 
provides some insight into the properties that could 
be expected of typical IACT GRB detection. In Ap- 
pendix |Bl we discuss how our results are affected if 
only the prompt phase of the burst emission is consid- 
ered, and the fading afterglow signal is ignored. 

4-1.2 Variation of model parameters 

In this section, we discuss the impact that variations 
in instrument properties and other general assump- 
tions could have on the GRB detection efficiency. This 
demonstrates the effect of variations from our baseline 
models discussed in the last section. A summary of re- 
sults is shown in Fig. [TO] 

The impact of VHE observation delay time due 
to GRB localization and telescope slew time, as dis- 
cussed in Section 12.3.31 is dependent on the assumed 
model for the GRB lightcurve at these energies. The 
upper left panel of Figure [TU] shows the overall impact 
of parameter Tdeiay on the detection efficiency, with 
other modeled parameters held constant. 

Next, in the upper right panel, we show how a 
higher or lower value of the telescope energy thresh- 
old than the ~ 25 GeV assumed in the baseline model 
would influence the detection efficiency. The effective 
area function is assumed here to have the same shape 
as presented in Figure [51 but with a shift in energy 
by a constant multiplicative factor. As discussed in 
the introduction, GRB observations are strongly af- 



GRB prospects for CTA 



13 




all GRBs 
CTA baseline, bandex 
CTA optimistic, bandex 



10000 
8000 
6000 
4000 
2000 


I 

>» 
o 

§j 0.5 
| 0.4 
™ 0.3 
1 0.2 
I 0.1 



all GRBs 
CTA baseline, fixed 
CTA optimistic, fixed 



4 

Redshift 



Fig. 9 The redshift distribution for detected GRBs in our model, for the bandex (on the left) model, and for the fixed model 
on the right. The upper panel on each side shows the number of detected GRBs when a baseline or optimistic (dotted) 
effective area function is assumed. The thin black line is the redshift distribution for all GRBs in the sample, which is 
created from the distribution shown in Figure [3] The lower panel is the fraction of detected GRBs in each bin. 



fected by spectral cutoffs due to EBL, and raising the 
telescope threshold energy reduces the redshift range 
over which GRBs are detectable. Detection efficiency 
is seen here to vary strongly with energy threshold, for 
both spectral extrapolation models. Note that setting 
the energy threshold here to 100 GeV is essentially the 
same as removing the LSTs from the telescope array 
(see Table ^ , since at these energies the effective area 
function is dominated by the MSTs. The large decline 
in detection efficiency with increasing energy thresh- 
old demonstrates the importance of having an LST 
array with low energy threshold to GRB detection, 
even though the LSTs may only contribute a fraction 
of the effective area of the total array at higher energy. 

The bottom- left quadrant of Fig. [TU1 addresses how 
altering the afterglow light curve index 7 in Eq. [5] af- 
fects results. As discussed in a]2.1.3l we have imple- 
mented a lightcurve in this work based on the T90 
time of a given GRB at Band peak energies, in which 
VHE emission is fiat for this period and then decays 
as £~ 3 / 2 . In such a model, 2/3 of the total VHE en- 
ergy emerges after the end of the T90 period, leading 
to a substantial afterglow flux that enables detection 
of GRBs after the lower energy emission has subsided. 
A faster or slower falloff of afterglow flux in time will 
change the optimal integration time for GRBs in our 
simulation, as well as the distribution in detection sig- 
nificance and therefore the detection efficiency. The 
effect is found to be relatively minor. 

Finally, we show in the bottom-right panel of the 
figure how altering the normalization of the back- 
ground rate changes detection efficiency. As discussed 
in ^2.3.21 the background rate assumed in this work 



is based on a rate of 6 photons per min over the 
VERITAS effective area (2 per min for the optimistic 
effective area), with extrapolation to lower energies 
achieved with a power law of index -2.7. This figure 
shows the effect of variations in this base rate. 

4-1-3 Spectral cutoffs and the impact of the EBL 



As mentioned in a!2.2.2l and the introduction, the EBL 
introduces a spectral cutoff in extragalactic gamma- 
ray observations that affects lower energies at higher 
redshift. In Table [3J we show how assuming differ- 
ent EBL models can change the overall detection ef- 
ficiency in our calculation. In general, the magnitude 
of UV/optical emission will determine the strength of 
the spectral cutoffs for GRBs at most redshifts. How- 
ever, at z £ 2 this UV emissivity of galaxies is highly 
uncertain by a factor of several (jGilmore et a 1 2009) 
(G09). The first two rows in Tableware the results for 
the GSPD11 EBL, which we assume elsewhere in this 
work. This EBL includes a UV-optical contribution 
that is nearly maximal in terms of the range allowed 
by high-redshift measured luminosity functions. The 
other rows in the table include models with less UV 
light; the fiducial model of G09 has a similar star- 
formation rate to GSPD11, but with more dust at- 
tenuation in high-redshift star-forming galaxies that 
reduces the UV emission. The low model in G09 has 
a smaller amount of star formation, in addition to the 



larger dust extinction. The model of lFranceschini et al 
(F08) predicts a similar UV flux to GSPD11 at 



low redshift, but a smaller amount at high redshift. 
Running our analysis with the low model of G09 and 
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Fig. 10 Effect on overall detection efficiency and detection rate of varying different assumptions about burst and instrument 
parameters. In each plot, the thick solid black and dashed blue lines are predictions from our bandex and fixed model, 
respectively, for the baseline effective area function. Thinner grey and cyan lines show the corresponding results using the 
optimistic effective area. The detection rates shown on the right-hand axes are based on an assumed instrument duty cycle 
of 0.1 for all telescopes, and an all sky trigger rate of 95 events/year; see i|4.3.1l Upper Left: The effect of varying the 
observation delay time (T^iay). Values are for the LSTs; the MST delay time at each instance is assumed to be the LST 
time plus 40 seconds. Upper Right: The effect on the overall detection efficiency and rate of changing the energy threshold 
of the baseline effective area function (Figure [6]) from its initial value of ~ 25 GeV. The optimistic area function differs from 
the baseline primarily due to a lower LST energy threshold (~ 10 GeV), and is not shown here. Lower Left: Variations in 
the assumed lightcurve slope index of the GRB afterglow (Eq. [2j, and the effect on overall detection efficiency and rate. 
Lower Right: The effect on the overall detection efficiency and rate of changing the background rate normalization by a 
multiplicative factor from its baseline value (parameter 'A' in Eq. 
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Table 3 Detection effic iencies found for a few different EBL 
models, including the lGilmore et all(|201lT ) (GSPDll) model 
used elsewhere in th is work, the fiducial and low models of 
Gilm ore et all (120091) fGO g). and the observational model of 
Franceschini et all |2008) (F08). For the F08 case, gamma- 
ray optical depths are only reported for z < 3 and we have 
used the z = 3 result for higher redshifts; this result is there- 
fore higher than it might be were the model extrapolated to 
higher redshift. The second column shows the effective area 
function assumed. 



the F08 model, we find detection rates of GRBs that 
are about 30 to 40 percent higher than the GSPDll 
case. However, the low model was disfavored in G09 
on the basis of IGM ionization data, in favor of the 
fiducial model, which only increases detection by 5 to 
15 percent over GSPDll. The F08 model only pro- 
vided gamma-ray opacities for z < 3, so opacities for 
high redshift GRBs may be artificially low in this case, 
and the detection efficiency inflated to some degree. 

Another caveat in our analysis is the possibility 
that GRBs typically have a spectral cutoff or turnover 
at some characteristic energy. Due to our lack of 
knowledge about the GeV-scale properties of GRBs 
beyond the few bright events that have been detected 
by Fermi and CGRO, it is difficult to explore this pos- 
sibility in detail. However, we can perform a simple 
test, and examine how our results change if a sharp 
spectral cutoff is assumed to exist at some characteris- 
tic observed energy. Obviously, if this energy is below 
the sensitivity region of CTA, then the detection effi- 
ciency must fall to zero, while if it is above the energy 
where the EBL has a strong impact on the spectra 
for the majority of detected GRB, then the effect on 
results will be minimal. In Figure [TTJ the impact of 
a universal step-function cutoff at a given observed 
energy is shown on the results for the total detection 
efficiency. 

This result suggests that our results will remain 
sound as long as emission continues unaffected to 
^ 100 GeV. In general, GRBs in the bandex model 
are found to be less strongly affected by this spectral 
cut, as detected bursts in this scenario include those 
with softer spectra than the universal spectral of -2 
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Fig. 11 The reduction in detection efficiency due to the in- 
troduction of a universal spectral cutoff at the indicated ob- 
served energy. The vertical axis of the plot is normalized to 
the detection efficiency found under our standard assump- 
tion that emission continues without a break to an energy 
of 1 TeV. The solid black line shows results for our bandex 
model, the dashed blue is for our fixed model. 



used in the fixed model. It is worth emphasizing that 
only GRBs with fluence greater than the median are 
generally detectable in our simulation (Fig. UHl , and 
so the existence of a fluence-dependent cutoff energy 
that affects only the fainter population of GRBs below 
100 GeV would have little effect on our results. 



4.2 Detection of Fermi-GBM bursts 



The GLAST Burst Monitor (GBM, iMeeean et al 
20091 )) on Fermi presents several unique challenges 
as a triggering instrument for ground-based follow- 
up. GB M is a potent sourc e of GRB alerts, about 
250 /yr (iPaciesas et alll2012l) . and if the Fermi mis- 
sion is extended to a 10-year period, ending in 2018 
or later, then there would be significant overlap with 
CTA operations and many alerts provided under op- 
timal viewing conditions for the array. 

Unfortunately, GBM is only able to provide ap- 
proximate coordinates for the GRB in real time, and 
the substantial uncertainties are typically similar to or 
larger than the field-of-view (FoV) of the LSTs, in con- 
trast to the arcminute localizations of the Swift-HAT. 
However, while the analysis of Section I4TT1 assumed a 
static observation centered on the target, one could 
also imagine strategies to compensate for the limited 
field-of-view, at the cost of exposure depth. After dis- 
cussing some specifics of our modeled observation of 
GBM alerts in the following subsection, we will ex- 



16 



R. C. Gilmore, et al. 



20 



Q0 
CD 
T3 



15 



O 
U 

H 10 

"cO 
o 



CO 

-t-> 

m 



—\ — i — i — r - 



- 1 — i — i — r - 



x*» 



■ v ,"",*,: v i x j 



_l I I l_ 



_l I I I I I l_ 



_J I I I I I I I l_ 



2-dimensional gaussian function, with total RMS an- 
gular uncertainty given by 



-8 -7 -6 -5 -4 -3 

Log Fluence [50 - 300 keV] (erg/cm 2 ) 

Fig. 12 Comparison of the fluences in the GBM burst popu- 
lation with the individual la statistical error on the position 
of each burst. 



plain one possible strategy for enhancing the DE of 
GBM bursts. 



4.2.1 Modeling GBM bursts 

Our properties for the population of GBM bursts are 
taken from a subset of 346 events from the upcoming 
GBM 2-year catalog (Paciesas et al. in prep). Data 
for these includes the statistical error on the position 
for each event, and the fluence in the 50 to 300 keV 
band. The relationship between the statistical error 
on the automated localization, produced in ground- 
processing and distributed within ~10 sec of the GRB 
trigger time, and the fluence in the 50 to 300 keV band 
is shown in Figure [T21 The 1 degree lower limit is due 
to the current grid size of the localization algorithm. 

For small statistical errors, a clear negative trend 
is seen between burst fluence and positional error. As 
the brightest GRBs are also the most easily detectable 
by CTA, this will work to our advantage in the detec- 
tion efficiency. GRB positional errors are also subject 
to considerable systematic uncertainty, which is in- 
dependent of the statistical error and shows no clear 
trend in brightness or other GRB properties at this 
time. We assume that 70 percent of the GRBs have a 
systematic error of 3.2 degrees, and 30 percent have a 
considerably higher value of 9.5 degrees (Connaughton 
et al. in prep). The total positional error then follows a 



Otot 



' a stat- 



(6) 



The fluence values available for this population 
cover a different energy range (50 to 300 keV) than 
the 20 keV to 2 MeV range covered by BATSE. Our 
model is based on fluences in the latter energy range. 
To compensate, we have calculated the ratios between 
the fluence over the full BATSE energy range and the 
50-300 keV range for the sample of BATSE GRBs. We 
then multiply the GBM fluences in the sample by the 
median of this collection of ratios, to produce a reason- 
able distribution of fluences over the BATSE energy 
range for the GBM population. The ratio found here 
is 2.14. 

A summary of fluence distributions in BATSE and 
GBM is shown in Fig.[T3] Recall that in Section PTTl we 
applied a multiplier of 0.501 to account for the differ- 
ences between the BATSE and Swift-BAT GRB pop- 
ulations. Our findings here suggest that GBM bursts 
are considerably brighter, on average, than those of 
BATSE, which are in turn brighter than Swift detec- 
tions (even the sub-set of Swift GRBs with known red- 
shifts). 

As GBM has substantially different energy cov- 
erage from Swift-BAT, and detects GRBs that are 
considerably brighter, one might expect the redshift 
distribution of GBM-detected GRBs to differ from 
the Swift-BAT population. Unfortunately, there are 
not enough known redshifts within the GBM popu- 
lation to do a comprehensive analysis of the differ- 
ences. This is largely due to the positional errors on 
GBM detections which make follow-up observations 
impossible without more accurate data from another 
experiment. We can, however, look at the handful of 
GBM bursts that are also listed with redshifts in the 
Swift catalogue. A plot of redshifts for these 23 GRBs 
is shown compared to the redshift distribution of all 
Swift GRBs in Fig. [TH GBM bursts are found to have 
a somewhat lower distribution of redshifts overall, and 
we have made a modification (dotted green line in the 
plot) to the fit for the Swift distribution that we will 
use in the analysis that follows. 

A final factor we must consider is an account of the 
camera sensitivity lost in off-axis CTA observations. 
We account for this by including a simple sensitivity 
factor that is a function of radius (r) from the FoV 
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Fig. 13 Integral distributions of fluence values for GRBs 
observed by BATSE and GBM, including GBM bursts be- 
fore and after the correction factor described in the text is 
applied. Solid red is the distribution of BATSE fluences over 
the full energy range of the experiment, 20 keV to 2 MeV. 
Dotted red are BATSE fluences in the 50 to 300 keV energy 
band, as determined by integrating over the Band function 
fits to the bursts. The green dotted line is the burst distribu- 
tion in the GBM sample, which is between 50 and 300 keV. 
The solid green line then shows the GBM sample after mul- 
tiplying by the median ratio (2.14) found for BATSE GRBs 
between these two energy ranges. Solid lines are therefore 
a direct comparison of the fluence distributions of the two 
instruments in our model. 




z 



Fig. 14 The integral distribution of redshifts in all Swift 
GRBs (blue) compared with the subset that also have GBM 
detections (green). Solid lines are the data for each, while 
the dotted blue and green lines are the fits used in the pre- 
vious sections and the following sections, respectively. 

center: 

S(r) = l;r<R fov -1.5 

S(r) = 0.3(R fov - 0.5 - r) + 0.7; 

R fov - 1.5 < r < R fov - 0.5 
S(r) = 0A(R fov - r) + 0.5; R fov - 0.5 < r < R fov 

(7) 



Here Rf ov denotes the maximum extent at which any 
observation is possible, and is equal to half the FoV 
value. This radial dependence gives full sensitivity up 
to Rfov — 1-5, 70 percent sensitivity at Rf ov — 0.5 
and 50 percent sensitivity at the edge of the camera. 
This factor is applied to both expected number counts 
from the GRB and the background rate, and is used 
for both LST and MST observations. We set Rf ov to 
4.25 degrees for the LSTs and 8 degrees for the MSTs. 

4-2.2 Static observations of GBM bursts 

As a first step, we calculate the detection efficiency of 
this population without applying any search mecha- 
nism in the GRB observation. This analysis, and that 
of the next section, will only utilize the fixed model, as 
spectral details for this burst population are not avail- 
able at the time of writing. Results from the previous 
sections indicate that detection efficiencies of GRBs 
are typically a factor of 1.5-2 times higher in the fixed 
model than in the bandex model; there is no reason 
to believe a similar relationship would not hold true 
as well here. 

The limiting factor in these observations is the field 
of view of the telescope. As most GRBs in the sam- 
ple have a minimum uncertainty of ~ 4 degrees from 
combined statistical and systematic effects, a 3 to 5 
degree diameter FoV is insufficient to catch more than 
a minority of events. Figure [T5l shows the current de- 
tection efficiencies calculated, and how these could be 
increased by future reductions in the amount of un- 
certainty affecting GBM burst positions as reported 
in real time. 

4-2.3 Scanning mode observations 

A possible solution to the FoV limitations of CTA is 
to attempt to rapidly scan over some portion of the 
GRB error box after the burst alert, rather than sim- 
ply observing the coordinates of the best GRB local- 
ization. We consider in this section the possibility of 
increasing the GRB detection efficiency using "scan- 
ning mode" observations, and address the question of 
how the search box should be chosen to optimize this 
rate. A larger search box increases the probability that 
the GRB will be observed, but at the cost of exposure 
time. We make the assumption here that the data from 
such a scanning observation over minutes or hours 
could be compared with a later determination of the 
actual burst position, and the significance of the GRB 
detection would then be computed in an after-the-fact 
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Fig. 15 Detection efficiency for static observations of GBM 
bursts, as a function of the uncertainty in burst position. 
Variations along the x-axis indicate the effect of multiply- 
ing the total positional uncertainty for each GRB, which is 
calculated as described in the previous section, by a given 
constant scaling factor. Curves show the detection efficiency 
for the CTA baseline (solid black) and CTA optimistic (dot- 
ted grey) effective area functions. The right-hand axis shows 
the detection rate for the telescope under a standard set of 
assumptions; see i)4.3l for details. 



analysis taking into account the photon counts and ex- 
pected background within the source PSF. We do not 
consider here the possibility that VHE emission could 
be identified in real time, i.e., for a "stop-on-target" 
type of scan. 

Our calculation of the detection efficiency follows 
that of i j4.ll above, with a few modifications. We as- 
sume that the time required for the LST array to 
search a region of 15 degrees radius is about 120 sec- 
onds. We therefore append the delay time used in the 
last section with another term that describes the de- 
lay between the commencement of the scan and when 
the telescope first passes over the GRB's true location. 
The total delay is then 



' delay 



= 60 + R N (0,120) 



(15°-R fov ) 2 



sec, 



(8) 



Rfov is one half the field-of-view for the LSTs, and 
i?Ar(0,120) is a random variable between and 120. 
Afterwards, it is assumed that the telescopes make 
many successive passes over the GRB, and that the in- 
tegrated exposure after Tdeiay is approximately equal 



to that of the standard calculation multiplied by a fac- 
tor (Rfov/Rsrch) 2 if the GRB is in the search box, zero 
otherwise. The received background is adjusted by the 
same factor. The MST array is assumed to search in 
the same pattern as the LSTs, and their total exposure 
is multiplied by an analogous factor. 

We consider several possible ways that one might 
determine the extent the scanned region, R sr ch- The 
simplest possibility is to use a constant for all GRBs. 
In a scan over possible values in steps of 0.5 degrees, 
we find that detection efficiency is maximized at R src h 
= 6 degrees for the baseline effective area, and 7.5 
degrees for the optimistic function. For the baseline 
effective area, the detection efficiency is increased to 
0.057 for R src h — 6 degrees, compared to 0.037 in 
the static case; an increase of over 50 percent. In the 
optimistic case, the detection efficiency can be more 
than doubled, from the static value of 0.07 to 0.16 
with an optimized scan box. These results show that 
a scanning strategy for GRB follow-up can raise the 
detection probability by a large factor, and that it 
is especially powerful for a telescope with a very low 
energy threshold. 



4.3 Total detection rate 

The second part of our calculation of the detection 
rate, summarized in Eq. [SJ is an estimation of the 
trigger rate from satellite instruments. This factor ac- 
counts for the sky coverage and duty cycle of the in- 
strument, and is the rate at which CTA can respond to 
and observe GRBs. All other factors influencing GRB 
detectability are incorporated into the detection effi- 
ciency parameter discussed above. Our calculations in 
the last section allow for observations at a maximum 
angle from zenith of 75 degrees, covering 37 percent of 
the sky. The calculations of detection efficiency include 
the effect of increasing energy threshold for observa- 
tions far from zenith. The duty cycle of Cherenkov 
telescopes is limited by the requirement that these 
telescopes operate mostly on clear, moonless nights, 
which has generally produced realistic values of about 
10 percent, and we expect this factor to remain valid 
for CTA. Operation of CTA during moonlight may in- 
crease the duty cycle to ~13 percent or more, but at 
the cost of a higher energy threshold. We therefore do 
not expect a significant change of our predictions for 
the overall detection rate. 
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The total trigger rate calculated from a given satel- 
lite alert rate 



TR = DC x SC x SR x BF. 



(9) 



Where the CTA duty cycle DC is - 0.1 and the sky 
coverage factor SC is 0.37, for a 75 degree radius area 
around zenith, as we have assumed throughout this 
work. The satellite rate SR is the number of GRB 
alerts produced by a given satellite detector per unit 
time. The bias factor (BF) includes the effect of any 
correlation or anti-correlation between the location 
of GRB satellite alerts and the sky area covered by 
CTA. This includes factors such as the bias present 
in Swift GRB alerts ( Gilmore et al 2010l) . which leads 
to GRB discovery preferentially in the anti-solar di- 
rection, working to the advantage of IACTs which are 
limited to nighttime operations. A factor of 1 indi- 
cates no departure from a random distribution of GRB 
alerts on the sky. We assume a factor of 1.4 for Swift- 
BAT GRBs and 1.0 for GBM. 

4.3.1 Swift-like GRBs 

The Swift satellite, launched in late 2004, has detected 
GRBs at a rate of about 95 /yr over its first 70 months 
of operation, and is expected to have an orbital life 
of >15 years ( Romano! |2010| ). If the science lifetime 



of Swift overlaps with that of CTA, then Swift will 
provide a constant source of well-localized GRB alerts. 
Using Eq. [§J we estimate a detection rate for Swift 
alerts of 



DR S 



DE x 4.92 



GRB 



If we assume the 'best-guess' instrument parameter 
of Tdeiay = 60 sec for the LSTs and use the detec- 
tion efficiencies from the first two rows of Table [3J we 
find detection rates of 0.37 and 0.57 yr -1 for the ban- 
dex and fixed models with the baseline effective area 
function, and 0.80 and 1.61 with the optimistic effec- 
tive area function. These correspond to timescales of 
32 and 21 months between GRB detections with the 
baseline area, and 15 and 7.5 months with the opti- 
mistic effective area functions. 

Another upcoming mission that could provide 
timely GRB localizations is t he SVOM satellite. As 
described in iGotz et all (|2009l) . SVOM will consist of 
an orbiting gamma-ray telescope covering an energy 
range similar to that of Swift that is expected to detect 
~70 GRBs/yr, as well as ground-based telescopes for 
follow-up observations. Because the ECLAIRs/CXG 



instrument on SVOM is intended to cover a similar 
energy range to Swift-BAT, the population of GRBs 
detected with this satellite should be similar to the 
Swift population, and therefore have a similar detec- 
tion efficiency to that which we have calculated here. 

4.3.2 Fermi-like GRBs 

GBM on Fermi has detected GRBs at a rate of about 
250 per year (Paciesas et al. in prep). Because these 
GRBs are detected at all points above the horizon, the 
anti-solar bias factor affecting Swift (and presumably 
SVOM) does not apply here. We can therefore write 
the detection rate for GBM GRBs as 
GRB 



DR G 



DE x 9.25 



yr 



meaning that GBM should provide around 10 alerts 
per year that can be investigated by CTA. While 
these GRBs are brighter on average than Swift GRBs, 
ground-based followup is hampered by the large un- 
certainty in burst location, which is generally several 
degrees. In M4.2.3I we showed that the detection effi- 
ciency of GBM alerts can be boosted by executing a 
rapid scan over some portion of the error box. The op- 
timal values found in our case for a fixed-model type 
of extrapolation, 0.057 and 0.161 for the baseline and 
optimistic effective areas, respectively, lead to typical 
detection timescales of 23 and 8 months. The values 
are similar to those found for the fixed model with 
Swift-BAT alerts. It is not possible to do a bandex- 
type analysis for these GRBs, as spectral information 
is not available at the time of writing. However, we can 
speculate that such a calculation would likely lead to 
detection efficiencies that are a factor of 1.5 to 2 lower 
than for the fixed model, which was the general finding 
for the BATSE GRB population. Therefore, we con- 
clude that Fermi-GBM and Swift-BAT could give rise 
to detection rates that are generally the same. How- 
ever, we point out that in the case of GBM alerts, the 
scanning mode necessarily means intermittent expo- 
sure on the source with therefore only a partial cover- 
age of the temporal emission of the detected GRBs. 

It is also possible that future improvements to real- 
time trigger analysis of GBM bursts could lead to bet- 
ter localization information. Significant improvements 
in this area could eliminate the need for scanning or 
other means to compensate for position uncertainty. 
As a simple test, we can take this possibility to an 
extreme and examine a case in which all positional 
uncertainty is removed from GBM alerts. In this case, 
the detection efficiency values are 0.19 and 0.45, which 
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Parameter 


Value (080916C) 


Value (BHLZ) 


Value (VA) 


Flux 


4.88 x lCT 3 


1.0 x 10~ 2 


8.8 x 10~ 5 


r 


-2.16 


-2.1 


-2.0 


T90 


66 s 


100 s 


50 s 


Redshift 


4.35 


0.5 


2.14 



Table 4 Parameters assumed in modeling the three GRBs 
described in jj5] These include the time-integrated flux 
normalization at 1 GeV (with units of GeV -1 cm -2 ), 
the spectral index in dN/dE, the T90 duration, and red- 
shift. Columns show parameters for GRB 080916C f€5TT)l. 
a "bright, hard, low-z" GRB ( tl5.2[) , and a "very average" 
GRB with parameters selected from the medians of the fixed 
sample used in this paper f <15.3l) . 



give detection rates of 1.8 and 4.2 GRBs yr _1 , ~ 3 
times greater than the optimal rates from our scan 
mode simulation. 



5 Results for specific GRBs 

It is useful to consider the spectrum that might be 
provided by an actual GRB detection. In this section 
we show sample spectra from a few different possible 
GRBs, which are modeled using the parameters sum- 
marized in Table 2) 



5.1 080916C 
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Fig. 16 A simulated realization of the detected spectrum 
from GRB080916C, assuming the parameters of Tableland 
an observation at 20 degrees from zenith. The top panel is 
for the baseline effective area function, and the bottom is 
for the optimistic function. The blue/grey points refer to 
the observed spectrum without any attenuation from the 
EBL. Black points are aft er applying the gamma-ray opac- 
ity of lGilmore et all (|201lh . Error bars shown only consider 
Poisson error in each bin. 



GRB 080916C was seen on September 16, 2008, by 
Fermi-LAT and GBM, soon after the beginning of sci- 
ence operations with the instrument. Th is GRB is no- 
table both for its high redshift (z = 4.35; Greiner et all 
2009) and its extremely high isotro pic-equivalent lum i- 
nosity, 8.8 x 10 54 ergs, or 4.9 M c 2 (|Abdo et alll2009df ). 
The finding of > 10 GeV emission from this burst 
can be used to set upper limits on the amount of UV 
light emitted from sta r - forming galax ies at high red- 
shift ()Abdo et ajl20ldlGilmorell201ll) . We can there- 
fore consider this GRB as an archetypical example 
of a bright, high- redshift GRB with a hard spectrum 
known to extend into the multi-GeV energy range. 

We model the high-energy emission from GRB 
080916C using the parameters shown in Table |H As 
in previous sections, an unbroken intrinsic power law 
extending to 1 TeV is assumed. The time-integrated 
flux and high energy spectrum are found using a time- 
weighted average of the spectra over GB M and LAT 
energy ranges as presented in Table 1 of Abdo et"al 
(|2009dl) . 

As in 21 the GeV lightcurve is assumed to decay as 
after the T90 period, with no spectral evolution. 



Following our analysis with these assumptions, we find 
that GRB 080916C could be detected at an angle from 
zenith as high as 39 degrees with the baseline effective 
area, or 58 degrees with the optimistic area function. 
In Figure 1161 we show the spectrum that could be ex- 
pected from an observation of the burst at a zenith 
angle of 20 degrees. In this figure, the spectrum is 
shown with a bin size of 0.1 dex, along with Poisson 
error bars for the total number of received photons 
(signal and noise) in each bin. For the baseline effec- 
tive area, we find a total of 682 signal photons received 
over the optimal integration timescale of 178 seconds. 
A signal is seen up to an energy of 90 GeV in each 
case, beyond which the signal to noise per bin is well 
below 1. For the optimistic effective area, the signal 
extends conclusively down to 10 GeV, and a total of 
17950 photons are detected over an optimal timescale 
of 562 sec. Note that the normalizations of the points 
in the two cases do not appear quite the same, be- 
cause integrated flux over T ^ s is being shown, and 
the timescales of integration are different. The effect 
of the EBL is easily seen in a comparison between the 
attenuated and unattenuated spectra, and the GRB 
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signal is discernible even at energies where an attenu- 
ation factor e~ T ~ 0.01 affects observations. 



5.2 Bright, hard, low redshift (BHLZ) 

We now consider the observed spectrum for a bright 
GRB observed at low redshift. For this "BHLZ" burst, 
we assume the parameters in the appropriate column 
of Table |H The very high GeV normalization of this 
GRB (about the 99th percentile in our bandex sam- 
ple) and its low redshift mean that it can be conclu- 
sively detected even at very large angles from zenith. 
In Fig. [T71 we show an observation of this GRB at an 
intermediate zenith angle (30°) and a very large angle 
(70°). In the second case, the energy threshold of the 
telescope is increased by a factor cos(70) -3 w 25, and 
therefore the observation is limited to energies above 
100 GeV. 



5.3 A "very average" (VA) GRB 

Finally, we repeat our analysis for a GRB that has a 
redshift, T90 duration, and high energy fluence chosen 
from the mean values of the fixed sample of GRBs (see 
Fig. [20] in Appendix [2]) ■ Properties are summarized 
the last column of Table @] Such a GRB is found to be 
only marginally detectable even occurring near zenith 
with a baseline telescope effective area. In Fig. [T51 we 
show the spectra from such a GRB. 



6 Conclusions 

In this work, we have attempted to make realistic pre- 
dictions of the GRB detectability with the CTA exper- 
iment, in particular by considering a reasonable range 
of possibilities for the GRB emission and the CTA re- 
sponse functions. The basic conclusions of this work 
can be summarized as follows: 

— CTA can be expected to conclusively detect one 
GRB every 20 to 30 months assuming a baseline 
effective area and background rate, or 1 to 2 GRBs 
per year with the optimistic instrument model. 

— Detected GRBs will be at a median redshift ~1, 
and a typical GRB detected with CTA will provide 
hundreds or thousands of signal events, mostly ap- 
pearing below 100 GeV. 



— Intrinsic spectral cutoffs will have little effect on 
our results provided that they are above an ob- 
served energy of 100 GeV (~ 200 GeV intrinsic). 

— The detection rate of GRBs is a strong function of 
the instrument energy threshold, and a somewhat 
weaker function of the typical response time. 

— Follow-up of GRB alerts from Fermi-GBM could 
benefit greatly from a scanning-type observation 
over the GBM error box, which can boost the de- 
tection rate by 1.5 to 2 compared to a static ob- 
servation. With such a strategy, the detection rate 
of Fermi-GBM bursts can be comparable to that 
of Swi/t-detected GRBs. 

— Bright GRBs will provide well-determined spectral 
information over at least a decade in energy, and 
will be a valuable source of information about the 
VHE emission mechanism and intervening cosmo- 
logical radiation fields. 

The detection rates we have determined are 
ro ughly in agreement w ith the independent estimate 
by Kakuwa et all ( 2011 ) that has been performed con- 
currently with our own work. Our findings are contin- 
gent on certain assumptions: namely that a satellite 
instrument (i.e., Swift, SVOM, or Fermi-GBM) will 
be available to provide burst alerts during CTA oper- 
ations at a rate similar to that seen in recent experi- 
ence, and that these alerts will be promptly transmit- 
ted and followed up with an instrument slew, when 
possible. These rates would approximately double in 
the case that two such satellites are available. 

Our predictions rely on a number of very uncer- 
tain assumptions about gamma-ray bursts that must 
be determined from limited data. These include the 
extrapolation of the spectrum to high energy, the typ- 
ical lightcurve of the high-energy component, and the 
amount of extragalactic background light which im- 
pedes observations of extragalactic sources in the GeV 
band. In the case of the first two, we have been guided 
by the observation by Fermi-LAT of £ 10 GeV pho- 
tons from 4 bright GRBs (080916C, 090510, 090902B, 
and 090926 A). A constant danger in this model is that 
these GRBs may not be representative of the popula- 
tion as a whole. If, in actuality, only a small fraction 
of GRBs have spectra that continue into the multi- 
GeV range, then our results here would overestimate 
the detection rate with CTA by a large factor. The 
aforementioned bright LAT-detected GRBs comprise 
less than 1 percent of the total observed by GBM and 
it is quite possible that spectral cutoffs routinely ex- 
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Fig. 17 A simulated observation of the BHLZ GRB, as described in the text and Table0] Spectra are for the burst observed 
at 30° from zenith (left) and 70° (right). Point types are as in Fig. 1161 Error bars shown only consider Poisson error in 
each bin. Note that the normalizations between the effective area functions are different because the optimal timescale of 
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ist between the GBM/BATSE energy ranges and the 
20-100 GeV band where most GRB photons would be 
detected (Figure [S]). Our bandex model incorporates 
some aspects of a cutoff for a significant fraction of 
GRBs; those bursts with (3 parameter — 3.0 have 



much less power in the GeV band than near the Band 
function peak, and are generally not detectable (see 
Fig. [21] in Appendix [A"| . Over one-third (36 percent) 
of GRBs in our model fall into this soft category. But 
given the limited energy range of BATSE, this num- 
ber is may not represent the full number of GRBs with 
spectral turnovers or cutoffs that are below the CTA 
energy coverage. 

Spectral turnovers or cutoffs could exist in GRB 
spectra due to internal absorption of gamma rays 
by source photons, or Klein-Nishina suppression of 
high-energy inverse- Compton emission that could be 
the ba s is for GeV-scale emission. As discussed in 
Barind (2006), internal absorption will lead to a spec- 
tral cutoff above an energy determined by the source 
bulk Lorentz factor (r). For r £ 1000 our results are 
likely unaffected. In general, only lower limits on f 
are available for GRBs; one possible exception being 
the bright LAT-detected GRB 090926 A, where the 
claim ed turnover in the GRB spectrum ( Breeeon et all 
20111 ) c an be interpreted as the effect of internal pair 
opacity ( Ackermann et alll2011 ) , and used to set limits 
on the bulk Lorentz factor: 200 ^ T ^ 700. Lower lim- 
its inferred from other bright LAT-detected GRBs are 
r > 900, 1200, and 1000 for GRBs 080916C, 090510 , 
and 090202B, respectively (jAckermann et a i 1201 ill . 
The distribution in r is generally unknown for dim- 
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mer GRBs, and it is possible that many bursts have 
factors in the hundreds, rather than thousands. How- 
ever, if it is the case that the brightest ~ 20 percent of 
GRBs tend to have unattenuated emission up to > 100 
GeV, then our results would not be strongly affected, 
as these are the events most likely to be detected in 
our model. 

It is unlikely that our results underestimate the 
CTA detection rate. An underestimate of brightness 
in the GeV band would entail an increase in the typ- 
ical MeV-GeV brightness ratio (Fig. [2]) above that 
seen for bright, hard GRBs with a LAT detection. 
The lack of LAT detections for most GRBs would 
seem to disfavor the possibility of large GeV-MeV ra- 
tios in fainter GRBs. Another possibility is that the 
EBL attenuation of the gamma-ray signal due to UV 
light for high-red shift GRBs is overes timated in the 
fiducial model of Gilmore et a 3 (1201 lh . We find that 
reducing the EBL flux to that of a minimal model, 
such as G09 Low in Table [3l could increase detection 
rates by as much as 30 to 40 percent, however disagree- 
ments between the G09 Low model and the bulk of the 
high-redshift data make this an unlikely possibility. Fi- 
nally we could have underestimated the performance 
of CTA itself, though the inclusion of our optimistic 
effective area model, with its sensitivity at energies as 
low at ~ 10 GeV and reduced background rate, is in- 
tended to be as hopeful as reasonably possible about 
the instrument capabilities. It is worth noting that 
CTA will produce a number of marginal detections, 
with 2 < a < 5. These are produced at roughly half 
the rate of bona fide (> 5a) detections in our models 
(Fig. [7]). These marginal detections would have higher 
significance than any of the GRB observations in th e 
VERITAS analysis fSection ETT1 lAcciari et all201ll) . 
where the highest significance quoted is 1.8. 

A basic assumption that we have made in this work 
is that GRBs to which the telescope is triggered, pre- 
sumably by a satellite instrument, will provide detec- 
tions at a much higher rate than serendipitous detec- 
tions, in which a GRB occurs in the field of view of 
the instrument during other observations. Simple ge- 
ometry shows that this is always the case. The sky 
coverage fraction of the LSTs, assuming a 4.25 de- 
gree field of view and, neglecting edge effects, is about 
0.14 percent. If we extrapolate from the GBM burst 
rate, and assume ~ 500 GRBs occurring per year over 
all sky, and include the duty cycle factor of 0.1, then 
we find 1 GRB inside the LST field of view during 
operations every 13 years. For the MSTs, with an 8 
degree field of view, we have one GRB every 4 years. 



Even with the enhanced detection efficiency that is 
possible with a delay time of zero (Fig. 1101 upper 
left panel), the detection rate for GRBs is still much 
less for the serendipitous case than for triggered op- 
erations. However, a serendipitous detection over the 
lifetime of CTA is possible, and would certainly be of 
great value in constraining the prompt VHE emission 
of the GRB. 

In H4.21 we discuss the prospects for detecting GRB 
afterglows from Fermi-GUM alerts. GBM bursts are 
found to be more than twice as bright on average 
than those from Swift, and a cursory look at the small 
number of redshifts available for this population sug- 
gests that they are typically slightly closer as well. 
The primary difficulty of detecting GBM bursts lies 
in the large positional uncertainty of the instrument; 
only about 10 percent of GRBs will fall in view of 
the LST for an observation at the center of the GBM 
error box, for a 4.25 degree field of view. We show in 
£ j4.2.3l how a scanning mode observation over the error 
box can improve the detection rate by £ 50 percent. 
With this change, we find that the detection rates of 
GBM bursts with CTA are nearly the same as Swift 
GRBs. If improvements to the GBM angular resolu- 
tion are possible before the onset of CTA, then our 
results could be enhanced by a significant factor, as 
described in Fig. [T5] We have considered the possibil- 
ity of improved localization with the use of a simple 
linear scaling of the positional error in this figure, as it 
is difficult to quantify to what extent such an error re- 
duction could take place over the next several years. It 
is worth emphasizing that because the brightest GRBs 
are generally those with the highest probability of de- 
tection, reductions in the errors for especially bright 
and/or hard GRBs will be the most advantageous in 
increasing the detection efficiency for GBM alerts. 

Our results show that GRB detection with CTA 
will rely heavily on the sensitivity achievable in the 20 
to 100 GeV band. As shown in Fig. [5J only about 20 
percent of the gamma rays found for typical detected 
GRB are at energies > 100 GeV; this is despite the 
large upturn of the effective area function in both of 
our assumed telescope models (Fig. [6]) at this energy. 
It should be emphasized that detection of GRBs with 
CTA is therefore heavily reliant on the performance of 
the LSTs. It should also be noted that without an LST 
component, we find that the detection rate for CTA is 
only marginally higher than for the VERITAS instru- 
ment (Table [5]) , and is less than one-third the rate for 
the complete CTA instrument. The importance of the 
LSTs is even greater if there exist intrinsic spectral 
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cutoffs in GRB spectra at energies between 20 and 
100 GeV. We note, however, that MSTs could play a 
significant role in providing high statistics data above 
100 GeV to perform time-resolved spectroscopy which 
could reveal much about the GRB physics. 

An exciting prospect is to use the spectrum 
of a GRB seen by CTA to constrain the EBL. 
While gamma-ray sources have been used in many 
cases to help cons train the EBL (see references in 
Gilmore et all 120111 ). these attempts have mainly fo- 
cused on relatively low redshift blazars. A high- 
statistics GRB detection by CTA at z > 1 could 
greatly improve our understanding of how the EBL 
evolves with redshift. Because the EBL impact is sig- 
nificant below 100 GeV at these distances, the LSTs 
will be crucial for such science. One great advantage 
of CTA is its ability to potentially detect simultane- 
ously both the attenuated and unabsorbed portions 
of a gamma-ray spectrum, which allows much more 
robust limits on the EBL than if only the attenu- 
ated spectrum is seen an d the intrinsic spectrum must 
be derived theoretically ( Raue and Mazin 2010l ). For 
GRBs at and above the median redshift of z — 1.7, 
an energy threshold of ^ 20 GeV will be needed to 
effectively capture the unattenuated slope. 
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A Other properties of detected GRBs 

In this Appendix, we examine in more detail how the pop- 
ulation of GRBs that pass our detection criteria compares 
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Timescale (s) T90 Duration (s) 

Fig. 19 Left: The distribution of integration timescales that maximize detection significance for detected GRBs, for the 
bandex (solid black) and fixed (broken blue) models. Right: Comparison of T90 for detected GRBs with the whole popu- 
lation, for 3 x 10 4 simulated GRBs. In the top panel, the thin line is the distribution of the full population, and the solid 
black and broken blue lines are the number of detected GRBs for the bandex and fixed models, respectively. The bottom 
panel shows the fraction of GRBs detected in each bin. 




Fluence @ 1 GeV (counts/GeV/cm2) Fluence @ 1 GeV (counts/GeV/cm2) 

Fig. 20 High energy fluence distribution of GRBs in our bandex (left) and fixed (right) models, together with the distribution 
for detected GRBs. Lower panels are the detected fraction of GRBs in a given bin. Solid black and dashed blue lines are 
distributions for baseline effective area, while dotted lines are the corresponding values for the optimistic area function. The 
thin solid lines in the top panels are the distribution for the full population. 



to the entire population of simulated Swift-like GRBs. This 
will give us some insight as to the properties that might be 
expected of a burst with a confident CTA detection. It will 
also be useful to look at how the assumption of different 
effective area functions can affect results. 

The distribution of integration timescales that maximize 
detection significance is shown in Figure [f9l The two spec- 
tral models produce similar results in this distribution. This 
result suggests that a integration timescale of 100 to 500 sec- 
onds after the commencement of ground-based observation 
will be favored for GRB detection in most cases, assuming 
a universal t -1 - 5 falloff in the afterglow lightcurve. A small 
subset of bright GRBs however are still visible against the 
background some hours after the event trigger (10 4 seconds 
in the longest timescale considered here), and nonzero re- 
sults are found for all bins in Fig. 1191 On the right hand 
side of this plot, we show how detection efficiency varies 



with GRB T90 duration, as determined by BATSE. Not 
surprisingly, longer bursts always have a better chance of 
being detected, but the majority of detected GRBs have 
T90 values from 30 to 100 sec, due to the scarcity of bursts 
with T90 > 100 sec. 

The differences between the fixed and bandex models 
become most apparent when we consider the distribution 
in high energy fluences predicted by each, as are shown in 
Figure 1201 In general, the bandex model has a much wider 
distribution in high energy fluence, because the beta pa- 
rameter introduces another degree of freedom into the ex- 
trapolation, and steep beta indices lead to a subset of the 
bursts in the sample having extremely low levels of high en- 
ergy emission. Conversely, the brightest bandex GRBs are 
brighter than the brightest bursts in the fixed model, as the 
latter are limited to a fluence ratio of 0.1 between ~1 MeV 
and ~1 GeV, while the corresponding ratio in the bandex 



GRB prospects for CTA 



27 




0.25 



0.20 



S 0.15 



0.10 



0.05 



0.00 




20 40 

Zenith Angle (deg) 



Fig. 21 The distribution of beta parameters from the ban- 
dex model. As in previous figures, the thin line is for the 
whole population, the thick solid line shows the GRBs de- 
tected using a baseline effective area, and the dotted line 
shows detections using an optimistic effective area function. 
The lower panel shows the fraction of detected GRBs in 
each bin. The rightmost bin in each panel designates GRBs 
that had f3 > —2 in the BATSE sample, and have been reset 
to -2 for this calculation. 



model can be as high as 1.0, with /3 = —2. This accounts 
for our somewhat unexpected result that while overall de- 
tection rates predicted by the bandex model are lower, de- 
tected GRBs in this model tend to be brighter than for the 
fixed model. Figure 1211 shows the distribution of /3 indices 
for detected GRBs in the bandex model. Only GRBs with 
fairly hard extrapolated spectra, /3 ^ — 2.5, are capable of 
being detected. 

Figure l2"2"l shows how the probability of detecting a GRB 
varies with the zenith angle zen at which it is observed. 
GRBs in our model are assumed to be observed at a single 
instantaneous point relative to zenith. While motion on the 
sky over the observation period T j, s will change 9 ze n over 
the course of longer integrations, the effect is small enough 
that we ignore it here. Detection efficiency is a weak function 
of 9zen out to ~40 degrees, where increasing energy thresh- 
old is compensated by increasing solid angles. At higher an- 
gles the detection efficiency declines more quickly. However, 
GRBs can in principle be detected out to angles as large as 
70 degrees, where the energy threshold is raised by a factor 
of 25 (Eq. [3]). These would have to be at low redshift, so 
as not to be completely obscured by EBL opacity combined 
with the elevated energy threshold of the telescope. 



Fig. 22 The fraction of GRBs detected as a function of 
zenith angle for the bandex model with baseline (solid) and 
optimistic (dotted) effective area functions. Results for the 
fixed model are qualitatively similar. 
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Fig. 23 Top: The integral distribution of photon counts 
arising from the prompt phase for detected GRBs in our 
bandex model. The solid line is for a baseline effective area 
function, while the dotted line is for the optimistic. Note 
that the y-axis intercept indicates the fraction of GRBs for 
which any photons are detected during the prompt phase; 
the majority of GRBs are detected purely on the basis of 
afterglow fluence. Results with the fixed model are quali- 
tatively similar. Bottom: The integral distribution of the 
fraction of high-energy fluence collected during the T90 pe- 
riod for detected GRBs; the remainder of the fluence being 
due to the burst afterglow. 



B Prompt phase observations 

GRB detection in our calculation is heavily reliant on emis- 
sion during the early afterglow phase. Only about 21 percent 
of GRBs in our sample have prompt emission (T90) phases 
longer than 60 sec, which we assume as a typical delay time 
for observations with the LSTs. The majority of GRBs are 
therefore completely inaccessible during the prompt phase 
for the standard assumption of a 60 sec time delay. As shown 
in Figure [T9l there is a definite bias toward longer duration 



GRBs in the detected portion of the population. Figure |2"31 
summarizes the amount of fluence in detected GRBs that 
arises from t < T90. About 57 percent of bursts detected 
with a baseline effective area have no prompt phase flu- 
ence, while only about 10 percent have more than half the 
detected fluence arising from emission during the prompt 
period. With an optimistic effective area function, the frac- 
tion of GRBs seen purely in the afterglow period is slightly 
higher. 



2S 



R. C. Gilmore, et al. 




Fig. 24 The integral distributions of detection significances (left) and photon counts for detected GRBs (right), in a scenario 
in which GRBs only emit during the prompt (T90) phase. Line types are as in Fig. 



We can also consider an extreme possibility in our de- 
tection efficiency calculation: one in which no high energy 
emission emerges after the prompt phase, or equivalently, 
the light curve index 7 in Eq.[3]is taken to +00. This is found 
to reduce detection efficiencies to about one-third their val- 
ues in the standard calculation: 0.027 and 0.040 for the ban- 
dex and fixed model with baseline effective area (0.057 and 
0.11 in the optimistic case). Fig. l24l shows the distribution of 
sigma values and counts for detected GRBs in such a case. 
These can be compared to those predicted in Fig. \7\ 



